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Abstract. We consider the problem of computing upper and lower bounds on the price of 
an European basket call option, given prices on other similar options. Although this problem 
is hard to solve exactly in the general case, we show that in some instances the upper and 
lower bounds can be computed via simple closed-form expressions, or linear programs. We 
also introduce an efficient linear programming relaxation of the general problem based on an 
integral transform interpretation of the call price function. We show that this relaxation is 
tight in some of the special cases examined before. 



Key words. Arbitrage, Linear Programming, Radon Transform, Basket Op- 
tions, Moment Problems. 

Acknowledgements. The authors are very grateful to two anonymous referees for numerous 
comments and corrections. We also benefited from discussions with Marco Avellaneda, Dimitris 
Bertsimas, Stephen Boyd, Rama Cont, Darrcll Duffic, Hans Follmer, Nicole El Karoui, Noured- 
dine El Karoui, Francesco Rapisarda and seminar participants at the groupe de recherche 
FIQAM meeting at the Ecole Polytechnique, the INFORMS 2003 annual meeting in Atlanta 
and the IMA workshop on Risk Management and Model Specification Issues in Finance. The 
work of the second author is supported in part by a National Science Foundation Career Award 
NSF-ECS-9983874. This paper first appeared as ERL Technical report UCB/ERL M02/33 in 
Nov. 2002 at U.C. Berkeley and on ArXiv as[math.OC/0302243]in Feb. 2003. 



Notation 

For two n-vectors x,y, x > y (resp. x < y) means Xi > yi (resp. Xi < y{), 
i = 1, . . . , n; x+ denotes the positive part of x, which is the vector with com- 
ponents max(zi,0). e is the n-vector with all components equal to one, and ei 
is the i-th unit vector of R". The set R™ denotes the set of n-vectors with 
non-negative components, and R™ + its interior. The cone of nonnegative mea- 
sures with support included in R™ is denoted by /C. For w G R m , K e R and 
g E R m+1 , the notation {g, (w,K)) denotes the scalar product g T w + g m+ iK, 
where g contains the first m elements of g. 
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1. Introduction 

1.1. Problem setup 

Let p G Rip, K Q g R+, w G R™ + and K { G R+, w l G R" , for i = 1, . . . , m. We 
consider the problem of computing upper (resp. lower) bounds on the price of a 
European basket call option with maturity T, strike Kq and weight vector Wq: 

max./min. 'E„(w r x — Kq)+ , 
subject to E 7T (wfx — K{) + — pi, i=l,...,m, 

over all probability distributions n G /C on the asset price vector x, consistent 
with a given set of observed prices pi of options on other baskets. Note that 
we implicitly assume that all the options have the same maturity, and that, 
without loss of generality, the risk-free interest rate is zero (we compare prices 
in the forward market). 

We seek non-parametric bounds, that is, we do not assume any specific model 
for the underlying asset prices; our sole assumption is the absence of a static 
or "buy-and-hold" arbitrage today (i.e. an arbitrage that only requires trading 
today and at the option maturity). The primary objective of these bounds is 
not to detect and exploit arbitrage opportunities in the basket option market, 
illiquidity issues are likely to make these opportunities hard to exploit. However, 
the data on basket prices (index options in equity markets or swaptions in fixed 
income) is very sparse and traders often rely on intuitive guesses to extrapolate 
the remaining points, using these prices to calibrate models and evaluate more 
complex derivatives. Our results aim to provide an efficient method to check the 
validity of these extrapolated prices where they are the most likely to create 
static arbitrage opportunities, i.e. very far in or out of the money. 

From a financial point of view, our approach can be seen as a one-period, 
non-parametric computation of the upper and lower hedging prices defined in 
|EKQ91| and |EKQ95| (see also |KS98| ). The necessary conditions we detail in 
scction0have been extensively used in the unidimensional case to infer informa- 
tion on the state-price density given option prices (see |BL78) or jLLOO] among 
others), we study here a multidimensional generalization. 

From an optimization point of view, problems such as the one above have 
received a significant amount of attention in various forms. First, we can think of 
problem P as a linear semi-infinite program, i.e. a linear program with a finite 
number of linear constraints on an infinite dimensional variable. We use this 
interpretation and the related duality results to compute closed-form solutions in 
some particular cases. Secondly, we can see QJ as a generalized moment problem. 
This approach was successfully used in dimension one by |BP02| . who solve 
the one dimensional problem completely and show that the multidimensional 
extension is NP-Hard. However, their relaxation algorithm requires the solution 
of a number of linear programs that is potentially exponential in n, the number 
of assets. This makes the method prohibitive for large-scale problems. Finally, as 
in |HS90| . one can think of as an integral transform inversion problem. This 
is the approach we adopt to design an efficient relaxation in the general case. 
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The contribution of this work is twofold. First, exploiting the necessary con- 
vexity of arbitrage free call prices 

C{w.K) = E n (w T x - K) + , 

where C{w,K) is the price of a basket call option with weights w and strike 
K, we detail a relaxation technique providing upper (resp. lower) bounds on the 
solution to . The resulting infinite dimensional linear program on the call price 
function can be solved exactly. Compared to the decomposition method proposed 
by |BPf)2| , this relaxation technique has the advantage of being polynomial-time 
in the number of assets and constraints. 

Secondly, in some particular cases, we provide exact solutions to that 
have a polynomial complexity in the number of assets and constraints. We also 
obtain expressions for the corresponding pricing measures, and use them later 
on to prove tightness of the linear programming relaxation in the general case. 
LWQ3j and |LW04| derived independently an equivalent upper bound in the 
same particular cases and exact lower bounds in dimension 2. 

Finally, using recent results by HLW04 who compute an optimal solution 
to Q in the particular case where one seeks an upper bound on the price of a 
basket option with positive weights given only single asset option prices: 

maximize 'E % (wqX — Kq)+ 

subject to E 7r (a; i - Kf) + =p[, i = 1, . . . ,m, j = 1, . . . , J\ 

we show that our linear programming relaxation is tight. 

Our paper is organized as follows. We begin in section H~^l bv a brief reminder 
of the fundamental duality between pricing measures and positive portfolios. In 
section |21 we detail a relaxation for problem using shape constraints on the 
call price C'(w,K) as a function of the weight vector w and the strike price 
K. Using linear programming duality arguments, we obtain in section closed- 
form formulas (or simple linear programs) giving upper and lower bounds in 
some particular cases. In section 0J using HLW04 we show that the linear 
programming relaxation derived in section [2] is tight in the particular case J5J) 
above. Finally, sectional provides some numerical examples. 



1.2. Semi-infinite programming duality 

We begin by detailing a key duality result linking the existence of a pricing 
measure (or state price density in Duf96 ) and the absence of an arbitrage 
portfolio. In the general case, we can write the upper bound problem as a semi- 
infinite program: 

p sup := sup / ip(x)Tr(x)dx subject to / (f)(x)Tr(x)dx — p, / ir(x)dx = 1, 
tteic Jr.™ Jr.™ Jr™ 

(3) 
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where 

ip(x) := (w£x- K )+, c/>i(x) := (wfx - Ki)+, i = 1, . 
We define the Lagrangian (on K. x R m+1 ): 



L(tt, A, Ao) = / ip(x)'K(x)dx+X [p— <p{x)iT{x)dx J +Ao(l— / n(x)dx), 
Jr™ \ Jr™ J Jr™ 

and, as in |HK93j . we can write the dual of as: 

d su P ._ inf . x r p + Ao . x T 4>(x) + A > ip(x) for every x G R" 

= inf : sup : A T p + ip(x) — \ T 4>{x). ^ ' 

A a;>0 

Both primal and dual problems have very intuitive financial interpretations. The 
primal problem looks for a pricing measure that maximizes the target option 
price while satisfying the pricing constraints imposed by the current market 
conditions. The dual problem looks for the least expensive portfolio of options 
and cash, A T 0(a;) + Ao, that dominates the option payoff ip(x). Of course, the 
dual problem above yields an upper bound on the upper bound. 
Similarly, the computation of the lower bound involves 

p lnl :— inf / ip(x)iT(x)dx subject to / (j}(x)n(x)dx — p, / n(x)dx = 1, 



R" 



whose dual is 



d inf := sup : X T p + A : \ T <t){x) + A < ip(x) for every x £ R" 

A ,A 

= sup : inf : \ T p + ij){x) — \ T (j)(x). 
A x >® 



(5) 



(6) 



Here, the dual problem provides a lower bound on the lower bound. 

General results on semi-infinite linear programs establish the equivalence be- 
tween the primal and dual formulations. We cite here a sufficient constraint 
qualification condition for perfect duality from |HK93| , which makes an assump- 
tion about the support of optimal distributions. (We focus now on the lower 
bound; a similar result holds for the upper bound problem.) 

Proposition 1. Assume that in problem the support of the asset price dis- 
tribution can be restricted to a given compact set B C R". Assume further that 
there exists a pair (Ao, A) G R n+1 such that: 

X T (j>(x) + A < ip{x) for all x G B. 

Then if d ini is finite, perfect duality holds, namely d lu! = p lnf . 
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Proof. See [HK93 . 



□ 



This constraint qualification condition trivially holds when 4>(x) and il>{x) 
are Call option payoffs hence we have d ml — p ia , provided that the support of 
distributions feasible for our problem can be restricted to some compact B C R™ . 
However, this is often not the case for the bounds detailed in section 13 and we 
will prove perfect duality directly whenever possible. 

2. Relaxation for the general case using an integral transform 

2.1. The Radon transform 

Let us come back to problem JU, for p G R™, K G R™, w G R" + , Wi G R™ , 
i = 1, ...,m and Kq > 0, we seek to compute upper and lower bounds on the 
price of a European call basket option with strike Kq and weight vector wq: 



with respect to all probability distributions tt G K, on the asset price vector x, 
consistent with a given set of m observed prices pi of options on other baskets 
and forward prices q% — E^Xj, that is, given 



Extending to basket options the results of |BL78j . we write, for some ir G JC: 
C{w,K) = Vt T {w t x~K)+ 

= I (w T x — K) + dir{x), 



we can think of C(w, K) as a particular integral transform of the measure tt and 
we can try to compute its inverse. If we assume that the measure tt is absolutely 
continuous with respect to the Lebesgue measure with density 7r(a;), then for 
almost all K we have: 



where 6(x) is the Dirac Delta function. This means that d 2 C(w, K)/dK 2 is the 
Radon transform of the measure tt (see Hcl99 for example). In this setting, the 
general pricing problem above can then be rewritten as the following infinite 
dimensional problem: 



Eniwox - K )+ 



E^ (wf x — Ki) + = pi, i=l,...,m and ~E 7T xj=qj 1 j = l,...,n. 





min./max. f(wo,Ko) 

subject to f(wj,Kj) =Pj, j = 1, . . . , m 
f(e i ,Q)=q l , i = l,...,n 

f{w,K)en c , 
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where ej is the Euclidean basis in R™ and IZc is the range of the (linear) integral 
transform 

C-.JC^TZc 

7r i— > C(w, K) = J Rra (w T .x — i4T) + d7r(a;). 

Thus, the problem of finding all possible arbitrage-free option prices becomes 
equivalent to that of characterizing the range of the Radon transform on the 



set of nonnegative measures /C. This has been done by HS90 in the context of 
production functions (which can be thought of as Put options). As in [HS90 , 
we denote by C°°{0, R"} the set of functions / such that for any k there is a 
polynomial Pk(x) of degree k such that: 

f(x) - P k (x) = o{\x\ k ) 7 as \x\ -> 0, x G R™ . 



Using Call-Put parity, we can directly derive from HS90, theorem 3.2] the fol- 
lowing result: 

Proposition 2. A function C(w,K), with w G R™ and K > belongs to IZc, 
i.e. it can be represented in the form 

C{w,K)= I {w T x - K) + dir{x), 



where n is a nonnegative measure on a compact o/R" , if and only if the following 
conditions hold. 

— C'(w,K) is convex and homogenous of degree one; 

— for every w 6 R™ + , we have 

lim C(w,K) = and lim ^ = -1; 

K^oo K^0+ OK 

— if we write = ^d/dxi, the function 

_ K ,fdC(w,K) 



FUr)= / . )K 

belongs to C°°{0, R™} and for some w G R" the inequalities: 

(-l) k+1 D ei ...D^F(Xw)>0 
hold for all positive integers k and A G R++ and all £i, in R" . 

Proof. See |HS90| . □ 



This result generalizes the necessary conditions for the absence of arbitrage used 
by |BL78j . jLLOOj or |BP02j in dimension one. 
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2.2. Linear programming relaxation 

The conditions above are not tractable in the general case but we can formulate 
a relaxation of the original program by replacing the last (moment) condition 
with weaker monotonicity and linearity conditions. We then get an upper bound 
on the upper bound (resp. a lower bound on the lower bound) solution by com- 
puting: 

sup/inf C(wq 1 Kq) 

subject to C(w, K) (jointly) convex in (w, K) 

C(w, K) homogeneous of degree 1 , . 

— 1 < dC(w, K)/dK < and C(w, K) nondecreasing in w 
C(Wi, 0) = wj q, i = 1, m 
C{wi, Ki) = pi, i = 1, m. 

where the variable is here C(w, K) G C (R" +1 — > R+). As we show below, this 
infinite program can be reduced to a finite LP. If we define p m +i — wj q and 
K m+ i = for i = 1, m and P2m+i — w^q with K2 m +i — 0, we can show the 
following result: 

Proposition 3. // the following finite LP: 
max./min. po 

subject to (gi,(wj,Kj) - (wi,Ki)) <pj -p^ i, j = 0, 2m+ 1 , , 

9ij>0, ~ 1 < 9i,n+i < 0, i = 0, ...,2m + 1, j = l,...,n 
(gi,(wi,Ki)) =pi, i = 0, ...,2m+ 1, 

in the variables po € R+ and gi £ R™ +1 for i = 0, 2m + 1, is strictly feasible 
and its optimal value is finite, the infinite program and its discretization 
10) have the same optimal value. Furthermore, an optimal point of TO can be 
constructed from the solution to 0). 

Proof. As in BV04 , we first notice that as a discretization of the infinite pro- 
gram the finite LP will compute a lower (or upper) bound on its optimal 
value. Let us now show that this bound is in fact equal to the optimal value of 
0. If we note z* — \pg, <?q T , . . . , g^ T ] the optimal solution to the LP problem 
above and if we define: 

C(w, K) = . max {p* + (g*, (w, K) - (w^K,))} , 

i— 0,...,2m-f 1 

where p* = pi for i = 1, . . . , 2m + 1. C(w, K) satisfies 

C(wi,Ki)=pi, i= 1, ...,2m + l, 

and, by construction, C(wq, Kq) attains the lower bound po computed in the 
finite LP. Also, C(w, K) is convex as the pointwise maximum of affine functions 
and is piecewise affine with gradient gi, which implies that it also satisfies the 
convexity and monotonicity conditions in 10, hence it is a feasible point of the 
infinite dimensional problem. This means that both problems have the same 
optimal value and C(w,K) is an optimal solution to the infinite dimensional 
program in 0. □ 
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3. Exact upper and lower bounds using LP duality 

In this section, we address the problem of computing exact bounds in some par- 
ticular cases using linear programming duality. We first consider the case when 
the observed data set corresponds to option and forward prices on each indi- 
vidual asset since in practice, observations always include the forward contract 
prices E^x* = i = 1, . . . , n (forward contracts exist whenever options do). The 
expressions derived in these simple particular cases will be used in the section 0] 
to show tightness of the upper bound relaxation in the general case. 

We first examine the problem of computing upper and lower bounds on 

E w (w T x - K )+, 

given the 2n constraints 

E w (xi -Ki) + =p h E 7r Xi = q i , i=l,...,n, (9) 

where Kq > and w,K,p,q are given vectors of R+_i_- We will assume that 
< p < q < p + K, which is a necessary and sufficient condition for the problem 
above to be feasible. We show sufficiency by constructing a discrete asset price 
distribution that matches these prices. Let us define marginal distributions TTi(xi) 
such that 

Jo with probability 7Tj(0) = 1 - 2^ 
Xl = \ S with Probability = ^ ( 10 ) 

and because 0<p<q<p + K, we know that 

0<«Z»<1. 

Ki 

Using Sklar's theorem, we can then construct an asset distribution n(x) with 
marginals ni(xi), i.e. matching the market prices of single asset options, hence 
the market is arbitrage free. From the form of the constraints in ©, we also 
observe that the constraints 0<p<q<p + K are necessary. 



3.1. Upper bound 

3.1.1. One forward and one option price constraint per asset Here, we apply 
the semi-infinite duality result to the upper bound problem described in @ to 
show the following result: 

Proposition 4. Let p, q G R" , Kq G R+, w G R"+ and Ki G R+ for i — 

1, . . . , n, with < pi < qi < pi + Ki. An upper bound on the optimal value of the 
problem: 

maximize E n (w T x — Ka)+ 
subject to E^x.;) = qi 

E w (xi - Ki)+ = pi, i = l,...,n, 
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is given by: 

(F up = max w T p + 7 u>,; whx(qi — pi, (3jKi) — /3jK , (11) 

0<j<n+l ^ — ' 
% 

with (3j := (qj —pj)/Kj £ [0,1], j — 1, ...,n and the convention [3q — 0, 
Pn+l = 1. 

Proof. In view of the general result Q), the dual problem can be expressed as 
d sup = inf sup X T p + [i T q + (w T x - K )+ - X T (x - K)+ - ^ T x, 

where, without loss of generality, we have included the constraint A + /i > w, 
in order to ensure that the inner supremum is finite. We introduce a partition 
of R™ as follows. To a given subset I of {1, . . . , n}, we associate a subset Di of 
RTf, defined by 

Dj = {x : Xi > Ki, i G I, < x. t < K l: i G I c } , 

where I c denotes the complement of I in {1, . . . , n}. For z G R", let Zj be the 
vector formed with the elements (^i)igj, in the ascending order of indices in I. 
We have for 

d sup = inf : max : sup : \ T p + ji T q 

\+H>w te{0,l},ZC{l,...,n} xGDj 

+t(w T x - Kq) - \J(xi - - [i T x 
= inf : max : X T p + ri T q + h(X, a, I.t), 

A+At>iu te{0,l}IC{l,...,n} 

where h(X, /i, /, t) is given by 

h{X, /1, /, t) :— sup : t(w T x — Kq) — Xj(xi — Kj) — /J 1 x 
= sup : (twin — fiic) T xic — tKo + XjKi 

0<x ic <K ic 

+ sup : (twi — iii — Xi) t xi 

(twic — fija^Kia — tK + (twi — fii) T Kj if Xj + flj > twi, 
-00 otherwise. 

We note that finiteness of h(X, fx, I, t) is guaranteed by A + fi > w and < 
t < 1. When these conditions hold, the maximum value of h(X, fi, I,t) over 
I C {1, . . . , n} is obtained when the complement I c is the full set, that is, when 
I is empty. We obtain 

max h(X, /1, /, t) = (tw — [l)\K — tKo. 

7C{l,...,n} 

Optimizing over i, we obtain 

max max h(X, u, I, t) = max ((—uS^K, (w — uS^K — Kn) . 
t£{0,l} ZC{l,...,n} vv ' + /+ > 
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This results in the following expression for d sup : 

d sup = inf X T p + ^ T q + max({-fi)lK, (w - (iffK - K ) 

X-\-fi>w 

= inf w T p + [i T {q - p) + max {(-^K, (w - nj^K - K ) , (13) 

which admits the following linear programming representation: 

d sup = inf :w T p + n T (q-p)+t t>v T K, v>0, v + ^i>0 

LA.t,V,Z 

t>z T K-K , z > 0, z + fi>w. 

The problem is feasible, and is thus equivalent to its dual. After some elimination 
of dual variables, the dual problem can be rewritten as: 

d sup = max w T p + w T y - (3K a : (1 - (3)K >q-p-y>0 
V ' P l3K>y>0. 

We remark that the above problem is feasible if and only if p < q < p + K. We 
thus recover the primal feasibility condition mentioned before. This condition 
ensures that the dual bound eP up is finite. The above further reduces to the 
one-dimensional problem: 

cf up = max : w T p + w, min(q l - p u (5Ki) - /3K . (14) 

— i 

The above problem is the maximization of a piecewise linear concave function 
of one variable, thus the maximum is attained at one of the break points 0j := 
(<7j ~ Pj)/Kj € [0, 1], j = 1, . . . , n, or for (3 = 0,1. This way, we can obtain a 
closed-form expression for the upper bound, namely 

(f up = max w T p + y~] Wi min(g i - p h (3jKi) - (3jK , 

0<j<n+l z — ' 
i 

with the convention /3o = 0, (3 n +i = 1, which is the desired result. □ 

Remark that, when the price of forwards is not given , the upper bound is readily 
obtained by setting the variable /i, which is the variable dual to the constraint 
E^a: = q, to zero in the expression (|13|l . We get the simple closed-form expression 

d sup = w T p + (w T K - K )+, (15) 

which can be obtained as a direct consequence of Jensen's inequality applied to 
the function x — > x+. We can check that the above bound satisfies some basic 
properties: it is convex in w and concave in p, q. Also, when w = (the i-th unit 
vector), and K = Ki, we obtain d sup = p i: while for K \ = (i.e. the options 
are in fact forward contracts), we obtain d sup = Pi(= 
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3.1.2. Two option price constraints per asset Using HLW04 , the result we 
just obtained can directly be extended to the (slightly) more general case where 
two option price constraints are given for each asset (but no forward price is 
specified) . We will use this technical result to show tightness of the LP relaxation 
in section 21 

We let p x ,p 2 e R£, K Q e R + , w G R!^ + and K 1 , K 2 6 R^ (with K 2 > K 1 ). 
We look for an upper bound on the optimal value of the problem: 

maximize E 7r (w T x — Ko)+ 

subject to E^Zj) = qi (16) 
E n ( Xi -Kf) + = p j i , i = l,...,n, j = 1,2. 

where q € R™ is here a variable. From J7J , to preclude arbitrage between options 
and forwards, we must impose: 

1 2 

£. := pj + Kl £iZ% <Qi<Pi+Ki :=qi, i = 1, . . . , n. 

As in |HLW04| . for each asset Xi, we can then form C^'(Kj), the largest de- 
creasing convex function such that C w (0) = p\ + K\ and C {l) (K\) = p j i for 
i = 1, . . . , n, j = 1,2. Then every maximal decreasing convex function C« (k) 
matching the market prices can be written in terms of qi as: 



From HLW04 , we know that for a given q the upper bound on l|16(l can be 
written: 

inf yVcW(^-jro) 

{A>0, \ T e=l\ ' Wi 
i— 1 

this means that an upper bound on the solution to (|16JI for all possible values 
of q can be found by solving: 

sup mf } Wl C Kl > —K - (pi + Kl - qi) —f 

{, <q t <q }{A>0, \T e=1} ^ \ Wj ) Kj 



i=l 



and because both C® and {Kl - (AlJT )) + are convex functions of A (with 
(pj + Kl — qt) < 0), this is also: 



inf sup ±w l C^(^K Q )-(pl + Kl- qi ) i -^^ :K0))+ 



{A>0, \ T e- 



The inner supremum is reached for qi = pj + Ki and we get the solution to l|16f) 
as: 



inf 



{A>0. A T e=l} 

z— 1 
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If qi = p\ + Ki, the corresponding measure in HLW04. places no weight on 
values of Xi smaller than K\, hence the problem with two options: 

maximize Et T (w t x — Kq) + 

subject to E^ixi - Kf) + =p{, i = 1, . . . ,n, j = 1, 2. 

is equivalent to the following problem (setting x = K 1 + y): 

maximize ~Ei 7T {w T y — (Kq — w T K 1 )) + 
subject to (z/i) = pj 

VAy* - (K? - K})) + =pl i = l,...,n. 

where one forward and one option price constraint are given per asset. 
3.2. Perfect duality: upper bound 

We first compute the optimal probability measures corresponding to the upper 
bound result with option and forward price constraints obtained in section l^.l.ll 
We can recover an optimal distribution, or a sequence of distributions which 
achieve the bound in the limit. This provides a direct proof of the fact that 
psup _ ^sup^ j e that the upper bound computed in Proposition 0] is tight. 
Without loss of generality, we assume e T w = 1. In l|14(l we obtained: 

d su P _ gup . w T p + w _ m i n / q . _ pi; - (3K , 

0</3<l i 

which can be rewritten (the min is taken elementwise): 

sup w T min{g — f3K e,p + (3 (K — K )}, 

0</3<l 

or again: 

sup inf w T {{e - t){q - (3K a e) + t{p + [3 {K t - K ))). 
o</3<i*e[o,i]>" 

Using LP duality we know that this is also equal to (with e T w = 1): 

inf sup f3(w T tK - K ) + w T (e - t)q + w T tp. 

t6[0,l] m n< / 3<l 

We express the above as 

inf w T (e ~t)q + w T tp+ (w T tK-K )^. 
te[o,i] m ' + 

This problem can be solved exactly as a finite linear program, and we obtain t* 
such that: 

d-P = W T (( e - t *)q + t*p) + (w T t* K - K ) . (17) 
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We recognize here the expression of the upper bound on the price of a basket, 
where we are only given the following option price constraints: 

E Tr (x i - Ki) + — pi, i=l,...,n, 

where K :— t*K and p := (e — t*)q + t*p. In other words, as in |HLW04| . this 
upper bound is equal to the bound we would obtain given only one (synthetic) 
call price pi per asset corresponding to a strike price Ki. 

Suppose first that (w T K < K ), in this case we have t* — e and d sup = w T p, 
and as in HLW04, p. 15, Case 3] we can construct a sequence of distributions 
converging to the optimal value. If however (w T K > K ) and p > 0, then using 
HLW04, p. 14, Case 2], we know that the bound is attained by a distribution with 
finite support. Finally, if (w T K > Kq) and pi = for some index i, HLW04 
p. 12, Case 1] shows that a similar result holds. 

Finally, using the result above, we know that this result is still valid when 
we replace the option and forward price constraints with two option price con- 
straints. 



3.3. Lower bound 

We can obtain a similar result for the lower bound problem. In this case however 
the solution is not in closed form and involves solving a (polynomial size) linear 
program. 

Proposition 5. Let p, q £ R" , Kq e R+ 7 w 6 R-++ and Ki £ R+ for i = 

1, . . . , n, with < p-i < Qi < pi + Ki. A lower bound on the optimal value of the 
problem: 

minimize E 7r (w T a: — Ko)+ 
subject to E^Xi) = cfe 

E^Xi - Ki)+ = Pi, i = l,...,n, 

can be computed by solving the following linear program: 
d ini = sup \ T p + fi T {q-K) + h 

X,fi,otQ,...,a n 

subject to A + fx < w 

h < a (w T K - K ) - (a w - fiftK, < a Q < 1 ( 18 ) 
Mi : h < ai(w T K - K a ) -^j^tiiaiWj - [ij) + Kj 
(A* + Hi)+/w l < cti < 1, 

which has 3n + 2 variables and 3n + 4 constraints. 

Proof. In the lower bound case, the dual problem is 

d int = sup inf \ T p + [i T q+ [w T x - K ) + - X T (x - K) + - fi T x, 

\+/J,<W x >° 

where we exploited the fact that the inner infimum is — oo unless A + /i < w. 
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Let us use the same notation as before. We have 

d mi — sup : min : inf : X T p + \j} q 

X+fj.Kw JC{l,...,n} xeD, 

+ {w T x - K Q ) + - \J( Xl - Ki) - fi T x 
= sup : min : X T p + /i T q + h(X, fj,, I), 

\+p<w /C{l,...,ra} 



where 



h(X, fj,, I) = inf ?/o - Xj (xi - K r ) - [i x : x e £>/, yo > w x - K , y a > 0. 
We have by linear programming duality 

h(X, fj,, I) — sup (aw — /i) T K — aKo — (awja — hjc^Kjc : awi — A/ — fj,i > 

< a < 1 

Thus 

d inf = sup X T p + (^(q -K)+ min /(A, fj,, I), 

X+p.<w JC{l,...,n} 

where 

f(X,fi,I) := sup a{w T K - K Q ) - (atB;« - /U/c)^^^, 

a(A, J u,J)<a<l 

and 

n (A J +/i J )+ 
a(A, jU, i) := max , 

iei Wi 

with the convention that a(A, fi, I) = when I is empty. 

Let 7 be a non-empty subset of {1, ... , n}. Let i G argmaxi S /(Ai + fii)+/wi. 
We observe that 

a(X,fJL,I) = a(X,fi,{i}), 

and 

/(A, M ,/)>/(A,M,W), 

which dramatically reduces the complexity of the minimization subproblcm: in- 
stead of computing the minimum over all 2™ sets I C {1, . . . , n} it is sufficient 
to pick I in the set of singletons of {1, . . . , n}, or I = 0. Therefore, the problem 
reads as a linear program 

d int = sup A T p + n T {q - K) + h 

\, n,a , ...,a:„ 

subject to X + < w 

h < a {w T K - K ) - (a w - n)\K, < a < 1 ( 19 ) 
Vi : h< a, L (w T K - K ) - Y,j^i( a t w j - Vj)+ K j 

(Xi + m)+/w l <cti<i, 

and can be solved efficiently, since it has 0(n) constraints and variables. □ 
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3.4- Perfect duality: lower bound without forwards 

Here we study a particular case of the lower bound problem discussed in sec- 
tion EH where we are not given information on forward prices. We can't prove 
perfect duality in the setting of section EH but prove it below in a more restric- 
tive case. Without information on the forward prices, we simply set the dual 
variable (i to zero in expression l|19fl to obtain: 

d ini = sup w T p£ + h : h < 0, h < ^(w.Ki - K ), 1 < i < n. f2 0l 
o<c<e v ' 

We note that d lni can now be expressed as the solution of a non-linear, convex 
optimization problem: 



d ini = sup w T p£, - max ^(K - w l K i ) + : < £ < e, (21) 

£ l<?<n 

or from its dual: 

n 

inf fa** ~ V <( K - w i R i)+)+ : » Te = !. v ^ °- ( 22 ) 



jinf _ 

v 

i=l 



We can reduce again this optimization problem to a line search over a scalar 
parameter, by elimination of the variable £. We obtain 

d ini = V" PiWi+ sup piiOj min(l, — — ^— — )-v. 

{i : KiWi>Ko} ~ {i : KiWi<Ko} 

The minimization above can be further reduced to a closed-form expression by 
noting that the piecewise-linear function (of v) involved has break points at 
ji = K Q — KiWi (for i such that 7, > 0) and 0. We obtain as before a lower 
bound on the minimization problem 

minimize 'E 7T (w T x — Kq)+ 

subject to E T (ii — Ki) + = pi, i — 1, . . . ,n, 

as: 



{i : KiWi>K } 



EKq — KjWj 
PiWi mil l, — — ) — Ko + WjKj 

\{l : KiWi<K } 

Now, the linear programming expression 11' I'D allows us to recover an optimal 
asset price distribution or a sequence of distributions that are optimal in the 
limit, as follows. Let v be an optimal vector for problem (|22|l . Let I be the 
set of indices i such that K > lUjifj. We note that i ^ I implies 2/j = 0. For 
simplicity we assume that I = {1, . . . , m}, where < m < n (the choice m = 
corresponding to empty Z) . 
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(23) 



First we examine the case when to = 0, that is, X is empty. In other words, 
mini w(i)Ki > Ko, and therefore d mf = p T w. For e > 0, we can define an asset 
price distribution 7r e such that 

e _1 p + K with probability 7r e = e 
with probability 7r = l — e 

which satisfies the pricing constraints and: 

E Tt [Kx - # ) + ] = u. T p + e(K - 

we recover the lower bound by taking the limit when e goes to zero. 
Next, we assume to > 1. Let a = (n — m)/m. For e such that 



e < a 1 min 1^(7^ 0), 

Kz<m 



we define the vector u{e) by 
i>i(e) 



— at if 1 < i < to, 
e otherwise. 



Since e is small enough, the vector z/(e) satisfies the constraints of problem (|22l) . 
We now define a distribution 7r e on the asset price vector x as follows. 

x = x e (i) with probability Vi(e), 

where 

'-Kjiij=i, 




x%i) 

otherwise. 

Note that Xj(i) is always well-defined, since Vj{e) > for every j. 

Let us check that the distribution 7r e of asset prices satisfies the constraints 
in Jy. For every j, 1 < j < n, we have 

E^fo - - ELi "itoO^W - 
= ^( e )(^(i)-^) + 
= Pj. 

We also check that with this choice of asset price distribution, the objective in 
ifTjl attains the lower bound d ml , when we let e — > 0. We have 

E^(w T x - if )+ - EjLi Vi[w T x*(S) - K ) + 

= EIU ^( £ )(E"=i wj^K*) ~ ^°)+ 

=Er=i^(e)K(^+^)-^o)+ 

= J2i=l( W iPi - v i( e )( K - WiKi)) + . 

Letting e — > 0, we obtain 

nm £ ^ E 7re (W T X - i^ )+ = E" - ^(-Kfj ~ U>t !£<)) + + E"=m+1 W iPi 

= E"=i( w i-Pi - ^(-^0 - w l K i ) + ) + 
= d int , 

as claimed. This shows that d lnf = p mf and that the lower bound computed in 
(|3.4(l is tight in the absence of constraints on forward prices. 
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4. Tightness of the integral transform based LP relaxation 

In this section, using results from HLW04 and section|21we show that the linear 
programming relaxation derived in Proposition |21 is tight, i.e. yields the exact 
solution to problem JJJ in the particular case considered in HLW04 where only 
market prices of single asset options are given for many strikes. 

Proposition 6. Let E R +7 K* S R+ for j — 1,...,J 4 , K G R+ and 
wq G R" + for i = 1, . . . , to. The problems 

maximize E„(wq x — l£o)+ 

subject to~Ejn{xi — Kj) + — p\, i = 1, . . . , n, j = 1, . . . , J x , 

and 

maximize pg 

subject to (g{, {uP^Kf) - (w^^fc)) < v\ - P\ hi = 0,...,n, k,l = 1,..., J 4 
ffi.i > 0, - 1 < < 0, k = i,...,n 

(9i, (wi,Ki)) =p h i = 0, . . . ,n, j = 1, . . . , J 1 

(25) 

where wf — e, for j = 1, . . . , J 1 arzrf J° = 2, /ia?;e i/ie same optimal value. 

Proof. We first focus on the particular case where only the forward price and 
one option price are given per asset, i.e. Kj = and J 1 = 2 for i — 1, . . . , n. 
As in section |3 we write qi = p| and in order for our problem to be feasible, we 
assume 0<p<q<p + K. 

Let us show that the LP in Ij25(l is feasible. Indeed, we can form a piecewise 
affine function that is feasible for Q by taking C(wo,Kq) = E w (wqX — Kq) + , 
where ir is the probability measure defined in l|10(l . By construction, this function 
corresponds to a feasible point of l|25(l and the variables gi in (|25|l are simply the 
sub-gradients of C(wq,Kq) at the data points. Finally, the LP in Q25JI is finite, 
since we always have < E 7T (wqX — Kq)+ < w$q and the feasible set of (|25|l is 
compact. This means that the optimum in l|25[) is attained. 

Suppose that the forward price information is ignored, i.e. assuming that 
to = 77, and wo G R™ . We note e*, the i-th unit vector. Without loss of generality, 
we set w^e — 1. Since the function C(wq,Kq) — w^p + (wqK — Kq) is a 
feasible point of the infinite LP in Q , if we call V LP the upper bound computed 
by the linear program (|25|l . we must have: 

V LP >w^p+(w^K -K ) + . 
Now, using the necessary conditions in (JZJ) and the convexity of 



(w T x - K) 
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in (w, K) we can write 

E Te (wlx - K Q ) + = (wlx - (w^K + (K - v$K))) + 

n 

< J2 w o,i®* c fa ~ fa + (Ko - WqK))) + 

i=l 
n 

i=l 

The conditions on the slope of the function C(w, K) imply 

n 

Wo,iC {e^K, + (K - u%K)) < w^p + (wqK - K ) + . 

i=l 

Hence, V LP < WqP + (wqK — Kq) and finally 

V LP =wlp+(wlK -K ) + , (26) 

where we recover the expression found in (|15|l . This means that the upper bound 
computed by the LP relaxation is tight in the particular case considered above. 

Now we turn to the case when forward price constraints E^Xi = q$ for 
i — 1, . . . , n, are included. As already observed in 13. 1.11 the function 

d sup (w , K ) = max : Wq p + V" w 0ji min(% - p h fyKi) - (3 3 K , 

0<j<n+l ^ — ' 
i 

is convex in (wq,Kq). Also, when wq = e,-, and Kq = Ki, we obtain d sup = p,-, 
while for Ki — 0, we obtain d aup — qi. This means that d sup (w, K) is a feasible 
point of the infinite program Q and hence V hp > d sup (wo, Kq). 
Since the finite LP H25|) is attained, at a point denoted by 

* r * *t *ti t 
z = [Po>9o ,---,9k \ . 

we can define the call price function 

d LP (u>, K) = . max { Pl + (g*, {w, K) - (w^K,))} , 

«=(),.. .,2n+l 

with po = Pq- We can compute the corresponding strike prices K = t*K and 
option prices p = (1 — t*)q + t*p, as in 13.21 By convexity of d JjP (w, K), we have 
d LP (ei, K) < §i for i = 1, . . . , n. From QZ®, we get that d LP {w 0l K ) = V LP < 
d sup {w 0l K ), hence finally d LP (w , Kq) = d sup (w , Kq). This shows that the LP 
relaxation of the upper bound is tight when the input is composed of forwards 
and one option price per asset. 

Now, from section r3.1.2l we know that the upper bound problem given two 
option prices per asset can be reduced to an upper bound problem given forwards 
and one option price per asset. This means that the LP relaxation of the upper 
bound problem is also tight in the case where two option prices are given for 
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each asset. Finally, we conclude by remarking that HLW04, Th. 4.1] show that 
the optimal upper bound involves at most two option prices per asset, hence 
the problem of finding an upper bound in the general problem l|24(l . i.e. given 
option prices for many strikes, reduces to an upper bound problem given only 
two options per asset, for which the LP relaxation is tight. □ 

The relaxation we obtain in section 2 treats the upper and lower bound prob- 
lems in a completely symmetric way. This does not at all reflect the inherent 
complexity of the exact problems. In the upper bound case, we are essentially 
lucky and we can find, in most cases, discrete distributions with compact sup- 
port matching the upper bounds, hence we very often get perfect duality (most 
notably in the special case discussed in this section). This is not true for the 
lower bound and we almost never get perfect duality in general. 

We can also remark that since the largest decreasing convex functions in- 
terpolating the market option prices are piecewise affine, the upper bound com- 
puted in HLWTmI Theorem 4.1] can also be formulated as the solution to a linear 
program. 

5. Numerical results 

In this section, we first show how a variant to the program in Proposition |3| can 
efficiently "clean" the market data on single asset option prices to remove the 
non convexities due to noise (asynchronous data, liquidity, etc). We then focus 
on an equity market example and compare the upper bound computed using 
the linear programming relaxation in Proposition and the bounds obtained in 
HLW04 on the data set described in |HLW041 Table 2] . We show that this same 
relaxation is not tight in the lower bound case by exhibiting an example where 
the lower bound computed in section f3. 31 is larger than the relaxation's result. 
Finally, we test our relaxation technique in the general case using a simplified 
interest rate data set taken from |BMf)l| and show how the lower and upper 
bounds perform in the general setting of computing bounds on the price of a 
basket given other basket prices. 

5.1. Basket bounds given single asset option prices 

5.1.1. Market data In Figure^ we show the market price on Mar. 17, 2004 of 
a certain number of basket call options on the Dow Jones index, with maturity 
Apr. 16, 2004, for various strikes (the underlying asset of this option is here the 
Dow Jones index divided by 100). We notice that the convexity requirement with 
respect to the strike price in Proposition's violated between strikes 95 and 103 
for example. Since prices in this plot are last quotes, this does not necessarily 
mean that an arbitrage is present. In practice, options with strikes far away from 
the forward price (here equal to 103) tend to be somewhat illiquid and their last 
emoted price can be an unreliable indication of the price at which they would 
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Basket Option Prices 

30 r 




80 85 90 95 100 105 110 115 120 

Strike Price 

Fig. 1. Dow Jones index call option prices on Mar. 17 2004, with maturity Apr. 16 2004, for 
various strikes. The vertical dotted line is at the forward price. 



trade today. However, since this option price data is used to calibrated pricing 
models, being able to efficiently detect such anomalies in the data is crucial. We 
can use the result of Proposition[3]to clean the input data on single asset option 
prices, from |LL00| we know that the relaxation in Proposition[3]is always exact 
in dimension one and we look for a decreasing, convex approximation of the 
market data. Let pi be the market prices of calls on an asset x with strikes Ki, 
for i — 1, . . . , m, we can find the closest (in the l\ norm sense) arbitrage free 
approximation of the market prices by solving: 

minimize Y%Li\Vi-Pi\ 

subject to (g i} (0, Kj - Ki)) < y 3 - y t , . 

<?;,!> 0, -l<. 9l , 2 <0 1 '> 

(gi,(l,Ki)) = y h i,j = l,...,m, 

which is a linear program in the variables gt E R 2 and j/j G R, for i — 1, . . . , m. 

5.1.2. Upper bound with many strikes Here we test empirically the result de- 
tailed in PropositionElusing the data set in [HLW04I Table 2]. This data set still 
contains minor non convexities (on the BA option prices for example) and we 
first clean it using the procedure described in l|27|l above. We then compare the 
upper bound on the price of a basket option with uniform weights 0.071 obtained 
using the relaxation in Proposition [3] with those detailed in HLW04, Table 3]. 
The results are detailed in Tabled As expected from the result in Proposition 
El the bounds match up to a small numerical error that is most likely due to 
differences in the data cleaning procedure. 
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Strikes DJX price LP relax. U.B. 



52 
56 
60 
64 
68 
70 
72 
76 
80 
84 
88 
90 
92 
94 
95 
96 
97 
98 
99 
100 
102 
103 
104 
105 
106 
107 



47.10 
43.10 
39.10 
35.10 
31.10 
29.10 
27.10 
23.10 
19.10 
15.20 
11.30 
9.40 
7.50 
5.80 
4.95 
4.15 
3.35 
2.73 
2.13 
1.60 
0.78 
0.50 
0.33 
0.15 
0.15 
0.15 



47.13 
43.14 
39.15 
35.16 
31.17 
29.18 
27.19 
23.21 
19.25 
15.37 
11.63 
9.80 
8.08 
6.49 
5.73 
5.01 
4.33 
3.71 
3.14 
2.64 
1.82 
1.48 
1.22 
1.00 
0.80 
0.64 



47.09 
43.10 
39.11 
35.11 
31.12 
29.13 
27.14 
23.15 
19.18 
15.24 
11.42 
9.61 
7.90 
6.32 
5.57 
4.85 
4.19 
3.58 
3.02 
2.53 
1.73 
1.42 
1.16 
0.95 
0.75 
0.59 



Table 1. This table displays for various strikes: the market price (DJX) of a basket option 
with uniform weights, the upper bou nd (LP re lax.) computed using the relaxation in j2land 
the upper bound (U.B.) computed in IHLW04I . 



Intuitively, tightness in this setting stems from the fact that the convex func- 
tion C(w,K) which is the optimal solution to the relaxation in Proposition [21 
coincides on each axis with the convex envelope of the single asset option prices 
in the data set (the functions (?W in |HLW04I p. 12]). Unfortunately, this does 
not happen in the more general case where the data consists of basket option 
prices. 



5.1.3. Lower bound We can also test the tightness of the various bounds ob- 
tained above on a simulated data set and compare these bounds with actual 
prices. The forward prices are given by q — {70,50,40,40,40}. We set K = 
{70, 50, 40, 40, 40} and p = {1.61, 1.43, 0.93, 0.70, 0.47}. Using the results of sec- 
tion |31 we get bounds on the price of a basket option with weight vector wq = 
{0.2,0.2,0.2,0.2,0.2} which are detailed in Table0 

We notice that the lower bound computed using l|19fl is tighter than that 
provided by the LP relaxation in © for at least one strike price. In this case, 
the LP relaxation is equal to the trivial lower bound given by the intrinsic value 
of the option. 
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Strike price 3.84 4.32 4.80 5.28 5.76 



Upper bound (relax.) 1.71 1.37 1.03 1.03 1.03 
Upper bound (from 1.71 1.37 1.03 1.03 1.03 
Lower bound (from © 0.96 0.48 0.09 0.00 0.00 
Lower bound (relax.) 0.96 0.48 0.00 0.00 0.00 

Table 2. Lower and upper bounds on the price of a basket given single asset forward prices 
and one option price per asset. We compare the lower bounds produced in J3] with those 
produced by the linear programming relaxation in J2] 



5.2. Basket bounds given basket option prices 



Here, we test the relaxation in Proposition[3]on a simplified data set taken from 
BM01 . Our objective here is to look for bounds on the price of a swaption given 
other swaption prices in a simplified setting. We refer the reader to |Reb98j and 
|d'A03j for further details. 

We compute market option pri ces by simulation. The underlying asset prices 
follow lognormal dynamics as in BS73|. The forwards are given here by F = 
{0.03, 0.04, 0.04, 0.05, 0.05} (these are forward interest rates), and the covariance 
matrix C is taken from BM01 pp. 301, 311]: 



C 



( 0.034 0.032 0.026 0.021 0.018 \ 
0.032 0.035 0.019 0.026 0.011 
0.026 0.019 0.024 0.010 0.019 
0.021 0.026 0.010 0.020 0.004 

V 0.018 0.011 0.019 0.004 0.017 / 



We look for upper and lower bounds on the price of a basket option (a yield 
curve option here) with weight vector wo = {.2, .1, .2, .1, .2} for various strikes, 
given the simulated price of other at the money basket options (swaptions) with 
weights vectors given by the lines of the matrix: 



W 



/ 1 0\ 
110 
1110 
11110 

\l 1 1 1 lj 



and the at the money price of a basket option with weights Wo (which is usually 
liquid). We are looking for bounds on basket options with weight wq and strike 
prices ±10% away from the forward. The resulting price bounds are plotted in 
Figure |2 in terms of implied volatility. We notice that the market volatility is 
actually equal (up to a small numerical error) to either the lower bound when 
the option is in the money or to the upper bound when the option is out of the 
money. 
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